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This  review  article  critically  outlines  and  discusses  the  main  issues  of  34  methods  which  have  been 
developed  and  validated  over  the  past  35  years  in  order  to  determine  with  an  acceptable  accuracy  and 
reliability  fundamental  parameters  of  solar  cells.  This  review  covers  methodologies  which  deal  with 
current-voltage  characteristic  {I-V)  analysis  either  theoretically  through  elaborated  models  and/or 
treated  graphically.  Methodologies  based  on  the  theoretical  analysis  of  the  I-V  characteristics  using 
the  one  or  two  diode  model  are  discussed.  The  investigation  on  the  l-V  characteristics  is  processed  via 
statistical  functions,  non-linear  regression  and  stochastic  models.  A  second  family  of  methods  to 
determine  the  solar  cell  electric  parameters  comprises  the  ones  which  deal  with  the  graphical  treatment 
and  analysis  of  the  I-V  characteristics  which  are  measured  at  different  environmental  conditions.  To  the 
third  family  belong  the  methods  which  use  a  mix  approach  of  theoretical  analysis  of  the  I-V 
characteristics  through  modeling  on  one  hand  and  the  graphical  analysis  of  their  experimental 
configuration,  on  the  other.  The  paper  comments  on  each  of  the  34  methods  and  provides  pros  and 
cons  for  the  determination  of  the  fundamental  electric  parameters  of  solar  cells. 
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1.  Introduction 

The  main  targets  in  the  energy  policy  and  the  research 
priorities  and  efforts  drive  towards  increasing  the  efficiency  of 
the  renewable  energy  systems  available,  and  to  the  design  of  cost 
effective  technologies. 

The  production  of  photovoltaic  energy  has  undergone  a  sustainable 
growth  rate,  with  maximum  points  in  the  three  oil  crisis,  when  the 
need  of  alternative  types  of  energy  sounds  prominent.  The  objective 
established  by  the  European  Union  is  that  until  2020  the  photovoltaic 
energy  represents  12%  of  the  total  electric  energy  production  [1], 
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The  determination  and  knowledge  of  the  influence  of  solar  cells 
parameters  upon  their  efficiency  become  crucial  for  the  optimiza¬ 
tion  of  fabrication  processes  and  for  the  scientific  research  [2],  The 
possibility  to  determine  the  important  parameters  of  the  solar 
cells  plays  a  major  role  in  evaluating  the  solar  cells,  in  controlling 
their  quality,  and  also  in  the  fabrication  and  power  performance  of 
reliable  solar  panels. 

The  most  important  parameters  of  solar  cells  can  be  deter¬ 
mined  by  using  the  current-voltage  (/-V)  characteristic  which  is 
shown  in  Fig.  1  and  by  analyzing  their  equivalent  circuit  [2j.  These 
parameters  are:  Iph  -  the  photogenerated  current,  lsc  -  the  short 
circuit  current,  Voc  -  the  open  circuit  voltage,  n  -  the  ideality  factor 
of  diode,  Rs  -  the  series  resistance,  Rsh  -  the  shunt  resistance, 
I0  -  the  reverse  saturation  current,  Pm  -  the  maximum  power,  Vm  - 
the  voltage  at  Pm,  Im  -  the  current  at  Pm,  FF  -  the  fill  factor  and 
>7  -  the  efficiency. 
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Fig.  2.  The  equivalent  circuit  for  a  solar  cell  -  the  one  diode  model. 


The  fill  factor  and  the  efficiency  of  a  solar  cell  can  be  calculated 


by  the  expressions 

PP  Pm  ImV  m 

~  IscV„c~  IscVoc 

(1) 

Pm 
^  ~  AIj 

(2) 

In  general  for  the  silicon  solar  cells  it  is  useful  to  take  into 
consideration  the  two  diode  model,  but  in  case  of  thin  film  cells 
(heterojunctions)  the  second  exponential  term  has  a  small  influ¬ 
ence,  which  can  thus  be  neglected  [5], 

2.  Methods  used 

The  measurement  of  the  I-V  characteristic  of  the  solar  cells  can 
be  done  in  the  laboratory  or  in  natural  sunlight  conditions.  In  the 
former  case  the  measurements  can  be  performed  under  dark  or 
illuminated  conditions,  in  the  latter  case  solar  simulators  being 
used  [6-25], 

Most  of  the  methods  used  to  determine  the  parameters  of  the 
solar  cells  use  the  one  diode  model,  as  the  interpretation  of  the 
equation  describing  the  mathematical  model  is  simpler. 
The  methods  to  determine  the  parameters  aforementioned  are 
briefly  presented  below. 

Method  1.  Chan  et  al.  developed  the  analytical  five  point 
method  [9],  This  method  uses  the  one  diode  model  for  solar  cells 
and  allows  the  determination  of  the  following  parameters:  Iph, 
n,  Rs  and  Rsh  (the  five  parameters  -  P)  under  illuminated  condi¬ 
tions  using  the  values  of  Isc,  Voc,  Im,  Vm,  the  slope  at  open-circuit 
point  Rso  and  the  slope  at  short-circuit  point  Rsho  which  are 
obtained  from  the  I-V  characteristic.  The  values  of  Rsho  given  by 
Eq.  (6)  and  Rso  by  Eq.  (7),  can  be  calculated  by  linear  fit  of  I-V 
characteristic  around  the  short  circuit  current  point  and  around 
the  open  circuit  voltage  point  respectively.  The  shunt  resistance  is 
calculated  using  the  following  expression: 

Rsh  =  Rsho  =  -  (  =  (6) 


where  A  represents  the  area  of  solar  cell  and  IT  the  solar  irradiance 
with  spectrum  AM  1.5  and  T=  25  C,  that  is  Standard  Test 
Conditions  (STC). 

The  solar  cell  electric  equivalent  circuits  and  the  corresponding 
mathematical  models  take  into  consideration  the  number  of 
diodes  in  the  circuit.  Fig.  2  presents  the  equivalent  circuit  of  the 
one  diode  model  which  is  mostly  utilized  to  characterize  the 
solar  cells. 

The  mathematical  model  which  describes  the  one  diode  model 
is  characterized  by  the  following  equation  [3]: 


f  —  Iph  —  lo  ( e  nVr  —  1 


V+IRs 

Rsh 


(3) 


Rs0  =  -(% 

dl2V  =  Va, 

The  ideality  factor  of  the  diode  can  be  calculated  by 
A 


Vt(B+C) 

where  A,  B  and  C  are  calculated  as  follows: 
A  =  Vm  +  RS0Im—V0c 


V , 


B=  In  (  hc-^-Im  )-  In 

Ksh 


V„ 


(7) 


(8) 

(9) 

(10) 


where  VT  represents  the  thermal  voltage,  VT=kT/ q,  k  denotes  the 
Boltzmann  constant,  T  is  the  temperature  of  solar  cell  and  q  is  the 
electronic  charge. 

In  case  the  electronic  conduction  mechanisms  within  a  solar 
cell  are  considered  separately,  this  can  be  modeled  by  a  power 
source  and  a  diode  in  parallel  for  each  mechanism  in  turn: 
diffusion,  generation-recombination  and  thermionic.  In  this  case 
the  above  mathematical  model  is  modified  by  adding  an  expo¬ 
nential  term  for  each  mechanism.  Thus,  the  two  diode  model, 
which  considers  the  diffuse  and  recombination  mechanisms,  is 
described  by  Eq.  (4)  [4],  while  Eq.  (5)  [4]  describes  the  three  diode 
model: 


1  —  Iph  Iqi 


e"iVT  -1 


V+IRs 

Rsh 


(4) 


where  d  and  r  indices  represent  the  diffusion  and  the  generated- 
recombination  mechanisms, 


1  —  Iph  fod 


V  +  IRs 

eW_i  )-/„ 


eW-i  -i  1 1 e+vr-1 


V+IRs 

Rsh 


(5) 


r _  Jm 

L  _  I  _)+ 

lsc  Rs* 

The  other  parameters  f0,  Rs  and  Iph 
ing  equations: 


.Voc 

nVi 


Rs  =  Rso -rlf^exp  I  - 


nV  7 


(11) 

are  obtained  using  the  follow- 
(12) 

(13) 


Iph  —  1st 


+  Io 


(14) 


Method  2.  Ishibashi  et  al.  proposed  a  new  method  for  the 
derivation  of  all  solar  cell  parameters.  The  parameters  are  obtained 
by  a  single  I-V  characteristic  measured  under  illumination  condi¬ 
tions,  using  the  one  diode  model  with  only  one  assumption, 
Eq.  (15),  which  is  true  for  the  most  types  of  solar  cells  [10], 


Voc-IscRs\  ;1 
nVT  J 


where  t  index  denotes  the  thermionic  mechanism. 


A  =  exp 


(15) 
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Using  the  equations  obtained  from  Eq.  (3)  for  the  coordinates  (0, 
Isc )  and  (Voc,  0),  and  the  assumption  Eq.  (15)  the  following 
expression  is  obtained: 


dV _  nVT 

~df  ~  lSc—I—[V—Rs(lsc—I)—nVT]/Rsh  +  5 


(16) 


,  _hc  Voc 
‘0  — 


,  y  =  exp — — 1 
y  r^sh  nVr 


P=V0C^-2nVT,  a=f§ 


b=l t 


dVj i  =  o 


(25) 

(26) 


Substituting  I=ISC  and  V=0  into  Eq.  (16)  the  relation  becomes 

=Rsh+Ks  =  Rsh  (17) 

where  Rs<tRsh.  Ishibashi  et  al.  introduced  two  temporary  para¬ 
meters  in  Eq.  (16),  Rso  and  n0  to  determine  the  series  resistance 
and  the  ideality  factor  of  diode.  These  two  terms  are  used  in  an 
iterative  process  [7],  By  plotting  -dV/dl  vs.  {/sc-/-[V-Fos(/sc-/)- 
n0Vr]/Rsll}-' ,  n  can  be  obtained  as  a  slope  and  Rs  as  ay-intercept. 
The  iterative  process  begins  with  Rso=  0  and  no=0  and  ends  when 
151  =  0,  where  5  is  given  by  the  following  formula: 

(Rs-Rso)(Isc-I)+(n-n0)VT 

V-Rso0sc-I)-noVT  1  ' 


Therefore,  the  other  parameters,  Ia  and  Iph  can  be  obtained 
using  the  following  expressions: 


lo  —  (  he  — 


V0c~hcRs 


Rsh 


exp 


__VV\ 

nVT) 


Iph  =  p-+Io  exp 


nVT) 


(19) 

(20) 


Method  3.  Chegaar  et  al.  developed  a  simple  method  to  extract 
the  following  parameters  of  the  solar  cell:  Rsh,  Rs,  Ia  and  n.  The 
parameters  are  obtained  under  illuminated  conditions  using  a 
single  I-V  characteristic,  the  one  diode  model,  an  auxiliary  func¬ 
tion  F(V)  and  a  fitting  routine  [2],  The  shunt  resistance  is  obtained 
using  Eq.  (6  ).  An  auxiliary  function  is  introduced  by  Chegaar  et  al. 
to  determine  the  other  three  parameters: 

F(V)  =  V-Va\n(Iph-I)  (21) 


The  assigned  particular  value  of  n  is  verified  by  comparison  with 
the  value  obtained  with  the  following  equation: 

_  (^m~^mRs){hc~^~h~^m~^m/Rsh)  t')l\ 

Vrilm-Vm/Rsh ) 

In  the  particular  case  when  a->  —  oo  Eqs.  (23)-(25)  become 
Rs=|b|-^,  Rsh^oc,  10JSC  (28) 

he  Y 

Method  5.  Easwarakhanthan  et  al.  developed  the  vertical 
optimization  method.  The  one  diode  model  was  used  for  solar 
cell.  This  method  allows  the  determination  of  the  five  parameters 
of  solar  cells  using  the  non-linear  least  squares  optimization 
algorithm  [12],  The  modified  Newton  method  with  the  Levenberg 
parameter  is  used  for  the  algorithm.  The  objective  function  S  in 
case  of  the  study  of  a  solar  cell  is 

S(P)=  2  (Ii-I(Vi,P))2  (29) 

i  =  1 

where  /,-,  V,  are  values  on  the  I-V  characteristic  whose  mathema¬ 
tical  model  is  given  by  Eq.  (3)  and  P  represents  the  set  of  the  five 
parameters.  The  goal  is  to  minimize  the  S  function. 

Method  6.  The  electric  conductance  optimization  method  was 
developed  by  Chegaar  et  al.  [13].  This  method  allows  the  determi¬ 
nation  of  the  five  parameters  of  a  solar  cell  using  the  non-linear 
least  squares  optimization  algorithm  considering  the  one  diode 
model  for  the  solar  cell.  The  method  is  based  on  the  modified 
Newton  method  with  the  Levenberg  parameter  used.  This  method 
considers  the  I-V  characteristic  and  the  slope  of  this  curve.  The 
objective  function  .S)  in  case  of  the  solar  cell  is 

S,(P)=  £  (G-C,(V„/,P))2  (30) 

i  =  1 


where  Va  represents  an  arbitrary  value  of  the  voltage.  Using  Eq.  (3) 
without  the  last  term  and  /  instead  of  V,  Eq.  (21)  becomes 

F(I)  =  aI+b\n(Iph-I)+c„  where  a  =  -Rs,  b  =  nVT  and  c  =  -nVTI0 

(22) 


The  series  resistance,  the  ideality  factor  of  the  diode  and  the 
reverse  saturation  current  can  be  determined  by  fitting  the  curve 
F(f)  vs.  I.  The  curve  is  obtained  by  calculating  F(f)  for  each  point  (V, 
/)  at  a  constant  temperature  and  for  a  value  of  Va.  Chegaar  et  al. 
take  into  consideration  integer  values  for  Va  between  1  and  5,  [2], 
The  photogenerated  current  is  considered  approximately  equal 
with  the  short  circuit  current. 

Method  4.  The  parameters  of  the  solar  cell  such  as  Rsh,  Rs,  Ia  and 
n  can  be  calculated  with  the  method  developed  by  Tivanov  et  al. 
[11],  Considering  the  mathematical  equation  for  the  one  diode 
model,  which  is  represented  by  Eq.  (3),  approximating  Iph  ss  he  and 
analyzing  the  I-V  characteristic  of  solar  cell  at  a  fixed  level  of 
illumination,  its  parameters  can  be  determined  by  the  following 
equations: 


(a“fa)2+B(a_fa)+G^ 


1/2 


+(a+b)+-y — 

he 


(23) 


Rsh  = 


nVT  nVr  i  , 
b+Rs  a  +  Rs^‘sc 


(24) 


where  C  is  the  conductance  calculated  with  Eq.  (31),  /,  and  Vt  are 
the  values  on  the  I-V  characteristic,  G,  represents  the  computed 
conductance  and  P  represents  the  set  of  the  five  parameters.  The 
goal  is  to  minimize  the  function. 


lph~h—(V +IRs)Gsh+nVTGsh 
nVT  +  Rs(Ipi,—I0—(V  +  IRs)Gsh  +  nVTGsh) 


where  Gsh  =  =- 
Rsh 

(31) 


Method  7.  Chegaar  et  al.  proposed  the  simple  conductance 
technique  based  on  Werner's  method  [12,13],  The  five  parameters 
can  be  determined  using  the  one  diode  model  and  Eq.  (3)  in  two 
steps.  In  the  first  step  the  shunt  conductance  Gsh  is  determined  by 
simple  linear  fitting  of  the  reverse  bias  characteristic  for  a  large 
negative  bias  voltage.  The  second  step  consists  of  rewriting  Eq.  (3) 
in  forward  conditions  for  V+RsI>kT  and  calculating  the  conduc¬ 
tance  G=dI/dV: 

c  _ 1  Rs  r  (32) 

lph-I  nVT  nVT  ( 


Using  the  plot  of  Eq.  (32),  the  ideality  factor  n  of  the  diode  can  be 
determined  by  the  y-intercept  and  the  series  resistance  Rs  by  the  slope 
respectively.  For  the  photogenerated  current,  it  is  highly  acceptable  to 
approximate  lph  k  1sc.  Taking  into  consideration  the  I-V  data  corrected 
using  the  effect  of  the  series  resistance,  Ia  can  be  determined  by  the 
standard  method  based  on  the  forward  I-V  data  [4[. 

Method  8.  Haouari-Merbaha  et  al.  developed  a  method  to 
determine  the  parameters  of  solar  cells  based  on  the  two  diode 
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model  [15],  The  I-V  characteristic  of  a  solar  cell  measured  under 
illumination  conditions  is  split  in  two  regions  by  the  resistive  line 
with  slope  V0JISC,  the  first  region  is  near  the  short  circuit  current 
and  the  second  region  is  near  the  open  circuit  voltage.  These  two 
regions  of  the  I-V  characteristic  are  fitted  as  follows:  the  first 
region  as  a  function  of  the  current  and  the  second  as  a  function  of 
the  voltage.  For  the  first  region  Eq.  (4)  can  be  modeled  as  Eq.  (33) 


and  for  the  second  as  Eq.  (34): 

l  =  lsc-GlhV-a3hs(I,V) 

(33) 

V=V0C-IR°-b3h°(I,V) 

(34) 

where  the  superscript  “s”  is  for  the  first  region,  the  superscript  “o” 
is  for  the  second  region,  h\I,V)  represents  the  curvature  of  the  I-V 
characteristic,  i  can  be  either  s  or  o  [14].  The  coefficients  (Isc,Gshs) 
can  be  obtained  by  fitting  the  first  region  near  short  circuit,  with 
Eq.  (33)  and  the  coefficients  ( V0C,RS° )  can  be  obtained  by  fitting 
the  second  region  near  open  circuit,  with  Eq.  (34).  The  set  of 
parameters  ( Iph ,  Gsh=I/RSh,  lod ,  Ior,  Rs )  can  be  determined  by  solving 
the  system  with  five  non-linear  equations  obtained  for:  Isc,  Gshs, 
Voc,  Rs°  and  the  curve  passes  through  the  boundary  point,  which  is 
around  the  maximum  power  point  [15],  The  ideality  factors  of  the 
diode  take  the  theoretical  values  nd=  1  and  nr= 2.  Solving  the 
system  starts  with  the  set  of  parameters  equal  with  zero  and  the 
iterative  process  is  done  when  the  parameters  converge. 

Method  9.  The  method  proposed  by  Ortiz-Conde  et  al.  to 
determine  the  parameters  of  a  solar  cell  is  based  on  the  Co¬ 
content  function  CC  [16].  The  one  diode  model  and  Eq.  (3)  are  used 
to  characterize  the  solar  cell.  The  I-V  characteristic  is  obtained 
under  illumination  conditions.  Using  the  Lambert  W  function  and 
the  Co-content  function  CC,  Eq.  (3)  can  be  rewritten  as 

CC(I,  V)  =  Cvl  V+Cn  (I -he)  +  Cn  vi  V(I-IX)  +  CV2  V2  +  C,2 ( I-Isc )2 

(35) 


The  equation  coefficients  Cvi,  Cn,  CVI  and  C/2  can  be  determined 
through  fitting  Eq.  (35)  using  the  CC  function.  The  CC  function  is 
numerically  calculated  from  the  measured  data  as  follows: 

CC(I,V)=  [V(I-Isc)dV  (36) 

Jo 


Provided  the  equation  coefficients  are  determined,  the  parameters 
of  the  solar  cell  can  be  calculated  by 


R 


1 

sh  ^  2Cv2 


(37) 


Rs  = 


y/I+I6Cv2Cl2-I 
4  CV2 


(38) 


given  by  Eq.  (15),  a  simplified  explicit  equation  can  be  obtained: 

'  =  J^Lamben  W  [nw  X)  eXP  („w)  *eXP„Vr  (Rs!sc  +  ] 


,  V  RshV 

Rs  sc  Rs(Rs  +  Rsh ) 


(42) 


It  can  be  observed  that  Eq.  (42)  is  an  explicit  equation  and  it  has 
only  three  parameters  n,  Rs  and  Rsh.  These  parameters  can  be 
determined  by  the  numerical  fitting  method.  The  initial  values  for 
these  three  parameters  can  be  obtained  as  follows:  the  shunt 
resistance  using  Eq.  (16)  and  using  the  linear  dependency  dV/dl  in 
function  of  (/sc+/-i//Rsh)_IkT/q,  the  series  resistance  and  the 
ideality  factor  n  of  diode  can  be  found  by  intercept  and  slope. 
The  saturation  current  and  the  photogenerated  current  can  be 
calculated  using  Eqs.  (19)  and  (20). 

Method  11.  Garrido-Alzar  developed  a  method  to  determine 
the  solar  cell  parameters  based  on  the  two  diode  model.  The 
method  consists  of  solving  the  equation  system  obtained  by  Eq.  (4) 
for  four  points  from  the  I-V  characteristic  ( V,-,  i  =  l,2,3,4  and 
nd=  1,  in  the  zero  order  approximation  for  n,  and  Rs  [18].  In  this 
step,  the  following  parameters  can  be  calculated:  Iph,  Iod,  Ior  and 
Rsh.  To  compute  nr  the  following  equation  is  used: 


Vt { In  [ Ior +IPh-Iod  (exp (%f- 1) ) -^]  -  In  Ior } 


The  series  resistance  is  obtained  using  the  following  equation: 


Is 


Iph  iod 


Vs+hRs 

Rsh 


(44) 


The  above  steps  must  be  repeated  until,  for  example,  the  value  of 
nr  is  obtained  with  the  desired  precision  [18],  Here,  nr^2  as  in 
method  8. 

Method  12.  Kiran  and  Inan  developed  a  method  to  determine 
the  parameters  of  the  solar  cell  based  on  the  one  diode  model 
under  illumination  conditions  [19],  They  obtained  an  approxima¬ 
tion  equation  from  Eq.  (3)  by  eliminating  the  saturation  current 
and  considering  the  shunt  resistance  as  infinite,  which  is  rendered 
as  Eq.  (45).  With  such  an  assumption,  it  is  clear  that  the  method 
provides  results  which  deviate  from  the  expected  values  when  the 
cell  is  aged.  This  has  a  value  in  the  field  of  study  of  the  PV  cell 
aging. 

V  =  V0C—IRs  +  mVT  ln|f^p/+exp;sc^v:^°c|  (45) 


There  are  only  three  parameters  Rs,  n  and  Isc  in  Eq.  (45).  These 
parameters  can  be  computed  by  fitting  the  experimental  data.  The 
photogenerated  current  can  be  obtained  as  Iph  ss  Isc.  The  saturation 
current  is  given  by  the  following  equation: 


Cv  i  (ymecv^ol-i)  +4  CnCV2 
4VjCy2 


(39) 


(Cy\  +  he)  (  \/I  +  16Cy2Cj2  +  1  ) 

Iph  = - V - 2 - L~2Cn  Cv2  (40) 

I—(V—IRs)/RSf,+Iph  .. 

0”exp((V-/Rs)/nVr)-l  ( 

Method  10.  Zhang  et  al.  developed  a  method  to  determine  the 
parameters  of  the  solar  cell  based  on  the  one  diode  model,  the 
single  I-V  characteristic  under  the  constant  illumination  level  and 
the  Lambert  W  function  [17],  The  implicit  Eq.  (3)  can  be  trans¬ 
formed  into  an  explicit  equation  with  five  unknown  parameters 
Iph<  I  O'  o,  Rs  and  Rsh  by  using  the  Lambert  W  function.  Using  Eq.  (3) 
in  the  particular  points  (0,  Isc )  and  (Voct  0)  and  the  assumption 


(46) 


Method  13.  Sellami  and  Bouai'cha  determined  the  five  para¬ 
meters  of  the  solar  cell  using  the  one  diode  model,  the  I-V 
characteristic  which  is  measured  under  illumination  conditions 
and  for  the  fitting  procedure  the  genetic  algorithm  (GA)  is  used 
[20].  The  GA  consists  of  some  preliminary  decisions:  solution 
encoding,  evaluation  function,  initial  population  generation,  selec¬ 
tion  criterion,  recombination/reproduction  and  termination 
criteria.  The  evaluation  function  used  in  GA,  derived  from  Eq.  (3)  is 

I(VuP)=Iph-I  o(e-^r-lJ— (47) 


The  cost  function  can  be  defined  in  the  case  of  solar  cells  as  the 
sum  of  the  squared  difference  between  the  experimental  and 
theoretical  current  values  given  by  Eq.  (29).  In  case  of  the  genetic 
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algorithm,  the  accuracy  criterion  being  given  by  the  cost  function, 
this  has  to  be  minimized.  The  cost  function  does  not  involve  the 
derivatives  and  it  is  determined  for  each  chromosome.  The  initial 
population  of  chromosomes  IPOP  can  be  calculated  as  follows  [12]: 


The  I-V  characteristic  is  measured  under  both  illumination  and 
dark  conditions.  In  order  to  reduce  the  number  of  fitting  para¬ 
meters  the  ideality  factor  of  diode  nd  takes  the  value  1.  The  series 
resistance  is  calculated  as  follows: 


IPOP  =  (hi-l0)  random  [Nipop,Npar\  +l„  (48) 

where  Npar  is  the  number  of  parameters,  5  in  this  case,  Nipop 
represents  the  number  of  chromosomes,  ft,  and  /„  are  the  highest 
and  the  lowest  values  of  the  five  parameters  respectively.  The 
selection  of  the  family  of  chromosomes  is  made  function  of  the 
cost.  The  genetic  algorithm  GA  is  stopped  when  the  value  of  the  S 
(P)  is  smaller  than  the  predefined  cost  minimum.  The  flow  chart 
described  by  Sellami  and  Bouai'cha  is  presented  in  Fig.  3.  The  same 
method  was  used  byjervase  et  al.  [21],  but  they  determined  the 
parameters  of  the  solar  cell  for  the  two  diode  model. 

Method  14.  Ye  et  al.  determined  the  five  or  seven  parameters  of 
the  solar  cell  using  the  one  or  two  diode  model,  the  I-V  char¬ 
acteristic  which  is  measured  under  illuminated  conditions  and  the 
particle  swarm  optimization  algorithm  (PSO)  [22],  The  same 
method  is  described  by  Qin  and  Kimball  [23],  The  PSO  algorithm 
allows  the  determination  of  the  solar  cell  parameters  without 
using  the  differential  equations.  This  algorithm  is  an  iterative  one, 
in  which  each  particle  flies  in  a  space  with  a  velocity,  which  is 
updated  as  a  linear  combination  by  the  position,  the  inertia 
weight,  social  iteration  and  self  cognition.  The  velocity  and 
position  for  each  particle  are  updated  as  follows: 

vi(k+\)  =  w(k)vi(k)+CrrVk)[p  best,(k)~x,(k)]  +c2r2lk)[g  best,(fc)— x,(k)] 

(49) 


Xj(f<+1)  =  Vi(k+1)+Xi(k)  (50) 

where  x,  represents  the  ith  particle  position  (each  particle  contains 
a  vector  with  the  solar  cell  parameters  values),  w  is  the  inertia 
weight,  Vi  is  the  velocity  of  ith  particle,  rj  and  r2  are  the  random 
numbers  between  0  and  1,  c,  and  c2  are  the  acceleration  constants, 
pbest  is  the  particle's  best  known  position,  gbest  represents  the 
best  known  position  found  by  all  particles  [22,23],  The  objective 
function  which  has  to  be  minimized  is 


F  = 


N 


k=  1 

N 


/f  =  |Wm| 


(51) 


where  Ip  represents  the  predicted  current,  Im  represents  the 
measured  current  and  ft  is  the  fitness  function  for  one  sampled 
point.  The  iteration  process  is  stopped  when  the  fitness  function 
reaches  minimum  conditions. 

Within  this  concept  of  solar  cell  analysis,  Sandrolini  et  al. 
determined  the  seven  parameters  of  a  solar  cell  using  the  two 
diode  model,  the  I-V  characteristic  which  is  measured  under 
illuminated  conditions  and  the  PSO  algorithm  and  cluster  analysis 
handled  together  [24], 

Method  15.  Stutenbaeumer  and  Mesfin  used  the  non-linear 
curve  fitting  to  determine  the  six  parameters  of  the  solar  cell  using 
the  simplified  two  diode  model  of  a  solar  cell  under  dark  condi¬ 
tions  [25]: 

/  =  /od(e*-l)+/0r(eV%-l)+^-  (52) 


The  simplified  two  diode  model  of  the  solar  cell  can  be  obtained 
for  small  values  of  Rs.  In  this  case,  the  approximation  Vt>IRs  holds. 


Fig.  3.  GA  flow  chart. 


Rs  =  - 


Vsc-Vol 


(53) 


where  Vsc  represents  the  necessary  voltage  in  dark  condition  to 
obtain  the  same  short  circuit  current  as  in  illumination  conditions. 

The  fitting  procedure  to  determine  the  parameters  uses  the 
non-linear  fitting  algorithm  developed  by  Levenberg-Marguardt 
which  is  implemented  in  Origin. 

Method  16.  Kaminski  et  al.  determined  the  following  para¬ 
meters  of  the  solar  cell  Iod,  Ion  nd,  nr,  Rs  and  Rsh  using  the  two  diode 
model  in  dark  conditions  [26,27],  The  I-V  characteristic  of  the 
solar  cell  in  dark  conditions  is  analyzed  on  two  regions:  the  first 
represents  the  higher  part  where  the  series  resistance  and  the 
diffusion  mechanism  are  important  and  the  second  represents  the 
lower  part  where  the  shunt  resistance  and  the  recombination 
mechanism  are  important.  The  first  part  is  used  to  determine  the 
following  parameters:  Iod,  nd  and  Rs.  The  series  resistance  and  the 
ideality  factor  of  the  diode  can  be  determined  using  the  linear 
regression  of  Eq„  (54). 


y  = 


1 

ndVT 


i-Rs-Fx)  > 


v=lnMI 
y  /-/,  • 


x  = 


V-Vt 

7=7T 


(54) 


where  ( V, )  is  a  point  on  the  first  part  of  the  solar  cell  I-V 
characteristic.  If  the  correlation  coefficient  is  not  good,  then  the 
value  of  I,  must  be  raised.  The  saturation  current,  can  be 
estimated  by  plotting  Inf  in  function  of  V-IRS.  The  second  part 
of  the  characteristic  curve  is  used  to  determine  the  following 
parameters:  Ior,  nr  and  Rsh.  The  linear  part  of  the  lower  part  of  I-V 
characteristic  is  fitted  by  the  least  square  method  using  Eq.  (55). 

/  =  /or[ew-ij  (55) 


The  shunt  resistance  is  calculated  then  by  Eq.  (56). 

f  lor  lor 

~  \dVj,  =  Is~ibVr~ihVT 


(56) 


Method  17.  Arcipiani  developed  the  generalized  area  method 
[28],  This  method  allows  the  determination  of  the  following  solar 
cell  parameters:  the  ideality  factor  of  diode,  the  series  and  shunt 
resistance.  The  method  consists  of  calculating  the  area  between 
the  OV  and  01  axes  and  the  I-V  characteristic  of  solar  cell,  Eq.  (57). 
The  one  diode  model  was  used  for  the  solar  cell,  using  the 
following  expression: 

A=  {Iph+Io){V0c-rIsc)r+Isc(\+gr)(^f~VTil\  +V0Cg(vTn-Vf\ 

(57) 


where  r=Rs  and  g=RSh .  Taking  into  consideration  the  following 
approximations  gr<  1,  l0<IPh  and  lsc&IPh  which  are  generally 
accepted  [29],  Eq.  (57)  becomes 


y  = 


2 wy+{v^VT"+iky{i)vrs"- y: 


IscVoc-A 

IscVoc 


(58) 


Eq.  (58)  contains  three  unknown  parameters  of  the  solar  cell,  r,  g 
and  n.  To  determine  these  parameters  Eq.  (58)  should  be  rewritten, 
taking  I-V  at  three  different  illumination  levels.  Then  the  resulting 
non-linear  system  of  equations  is  solved. 

Method  18.  The  method  developed  by  Jia  and  Anderson  allows 
the  determination  of  the  series  resistance  and  the  ideality  factor  of 
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diode  [30],  In  this  method  n  is  considered  variable  on  the  I-V 
characteristic,  thus  n  =  l  at  Voc  and  n  =  2  at  Isc.  The  one  diode 
model  is  used  and  the  solar  cell  is  illuminated.  The  series 
resistance  can  be  calculated  with  Eq.  (59)  which  is  obtained  using 
the  following  approximations:  Iph  ss  lsc,  Rsh  is  infinite  and  ISC>I0. 

„  Vm{-t/VT)(lsc-Im){Voc  +  VT  ln[l -{Im/hc)]}-Im 
s  lm(-l/VT)(Isc-lm){Voc+VT\n[\-(Im/lsc)]}+Im 


The  value  of  the  ideality  factor  of  the  diode  at  maximum  power 
point  can  be  calculated  as  follows: 

n  =  (Vm+ImRs)/(v0C+VT\n^c^\  (60) 


Method  19.  The  integral  method  was  developed  by  Ortiz- 
Conde  et  al.  [31]  and  Kaminski  et  al.  improved  this  method  [26], 
The  series  resistance  and  the  ideality  factor  of  the  diode  can  be 
obtained  under  dark  conditions  through  the  linear  regression  of 
Eq.  (62).  This  equation  is  obtained  from  Eq.  (61)  through  integra¬ 
tion.  The  one  diode  model  was  used  for  the  solar  cell  and  the  shunt 
resistance  is  considered  infinite. 

/  =  /0(e^r-l)  (61) 


(v 

y  =  nVT+Rsx,  x  =  '-±±,  y  =  TAK  /  IdV  (62) 

Jvx 

where  (\7? ,/j)  is  a  point  of  the  solar  cell  I-V  characteristic  and  /  is 
the  dark  current.  The  integral  can  be  solved  using  the  means  of 
trapezoidal  method. 

Method  20.  El-Adawil  and  Al-Nuaim  proposed  a  method  to 
determine  the  series  and  shut  resistances  using  the  one  diode 
model  [32],  The  series  resistance  can  be  estimated  by  using  Eq.  (3), 
with  the  approximations  shown  in  Eq.  (63)  and  by  choosing  two 
points  around  the  knee  of  I-V  characteristic  of  solar  cell  (V),  /,■), 
i=l,2. 


I0Rshe^>V^  §£</sc+/o-/i,i=l,2,  hc*>I0,  Isc1*  I, 


ph 


(63) 


Rs  = 


nVT 

h-h 


In 


Uc-h 

Isc-h 


V2-V, 

h~h 


(64) 


The  shunt  resistance  can  be  obtained  by  the  following  equa¬ 
tion: 

Rsh+Rs  _V\  ~Rs  .cir, 

nV,h 


Thongprona  et  al.  had  used  the  same  method  before,  but  for 
both  dark  and  illuminated  conditions  [33]. 

Method  21.  The  ideality  factor  of  the  diode  and  the  reverse 
current  can  be  determined  by  using  the  semi-logarithmic  fsc-Voc 
characteristic,  as  in  the  equation  below.  The  short  circuit  current 
and  the  open  circuit  voltage  must  be  determined  for  different 
levels  of  illumination.  In  this  method,  the  one  diode  model  is  used 
with  the  following  approximations: 

Voc 

Rsh^oo,  Isc  *  Iph,  eW*  s»l  (66) 


Under  these  conditions,  Eq.  (3)  becomes 

lti  Isc  =  In  I0  4-  V oc  (67) 


The  ideality  factor  of  the  diode  can  be  calculated  from  the  slope 
of  the  semi-logarithmic  characteristic  and  the  saturation  current  is 
calculated  using  the  intercept  of  the  characteristic  on  Y-axis. 

Priyanka  et  al.  improved  this  method  by  taking  into  considera¬ 
tion  the  value  of  the  shunt  resistance  [34  .  Eq.  (67)  becomes 


Method  22.  The  method  developed  by  Signal  allows  the 
determination  of  the  series  resistance  and  the  ideality  factor  of 
diode  using  the  one  diode  model  for  one  level  of  illumination  [35], 
The  shunt  resistance  for  this  method  is  considered  infinite.  The 
author  developed  analytical  relations  for  Pm,  Im,  Vm  and  FF  in 
function  of  n  and  Rs.  The  values  for  n  were  chosen  1  or  2.  The  series 
resistance  was  considered  ideal  (Rs=0)  or  real  {RSV=  0).  For  two 
values  of  n  and  R5i  A  FF  and  A  (Ym/Yoc)  were  calculated.  If  n  is  not 
equal  the  correct  value  of  n,  the  value  of  the  series  resistance 
obtained  from  A  FF  differs  from  that  obtained  from  A  (Vm/Voc).  The 
series  resistance  and  the  ideality  factor  of  the  diode  are  found  by 
varying  the  value  of  n  until  the  two  values  of  Rs  obtained  from  A  FF 
and  A(Ym/Yoc)  are  almost  equal. 

Method  23.  The  series  and  shunt  resistance  can  be  determined 
using  the  one  diode  model  and  the  I-V  characteristic  of  a  solar  cell 
in  the  third  and  fourth  quadrants  under  illumination  conditions. 
This  method  was  developed  by  Priyanka  et  al.  [34].  The  series 
resistance  can  be  calculated  using  the  following  equation: 


Rs  =  nV'\n 


nrR-(Vr+Vf+IfR)I 
Io  R 


Yl 

If 


(69) 


where  R=Rs+Rsh,  I/,  V/  represent  I,  V  in  the  fourth  quadrant  and  Ir, 
Vr  represent  /,  V  in  the  third  quadrant.  The  shunt  resistance  can  be 
determined  using  the  Eq.  (6)  and  by  making  the  approximation 
RS<R,  R  is  determined.  The  ideality  factor  of  the  diode  and  the 
saturation  current  must  be  previously  determined  using 
method  21. 

The  series  resistance  can  be  obtained  more  accurately  using 
Eq.  (70)  by  iteration.  The  first  value  of  the  series  resistance  in 
iterative  process  is  given  from  Eq.  (69). 


Rs  =  nV'ln 


IrR-(Vr+Vf+IfRy 

Io(R—Rs) 


h. 

If 


(70) 


The  values  for  /r  and  If  should  be  chosen  so  that 

IrR-(Vr  +  Vf+IfR)>  0  (71) 


Method  24.  Bowden  and  Rohatgi  developed  a  method  which 
allows  the  determination  of  the  two  parameters  of  the  solar  cell,  n 
and  Rs  [36],  In  this  method,  two  I-V  characteristics  are  used,  one 
measured  for  one  sun  illumination  and  the  other  measured  for 
0.1  sun  illumination  or  in  the  ideal  case,  at  such  an  illumination 
level  that  the  short  circuit  current  of  the  solar  cell  equals  the 
difference  hc.i—Im.t ■  The  series  resistance  can  be  calculated  using 


Rs  = 


V  oc.0,1  —  VmA 
he ,  1— fsc.0,1 


(72) 


where  the  pair  (Ymi,  ISc,t )  is  obtained  from  the  I-V  characteristic 
measured  for  one  sun  illumination  and  the  pair  (Vocai.  4c,o,i)  is 
obtained  from  the  I-V  characteristic  measured  for  0.1  sun  illumi¬ 
nation.  The  ideality  factor  of  diode  can  be  calculated  using 


v  oc,\—  ^oc,0,l  1 

ln(/Sc,l//sc,0,l)V/T 


(73) 


Method  25.  The  series  resistance  can  be  determined  by  the  two 
characteristics  method  which  was  suggested  by  Swanson  [37], 
described  by  Wolf  and  Rauschenbach  [38]  and  reiterated  by 
Bashahu  and  Habyarimana  [39],  In  this  method  the  I-V  character¬ 
istic  of  a  solar  cell  is  measured  at  the  same  temperature  for  two 
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levels  of  illumination.  The  series  resistance  can  be  calculated  using 
the  following  equation: 


AV  =  V2-V, 
s  A/sc  h-I2 


(74) 


The  series  resistance  can  be  determined  as  follows: 

R  Vm  1 

Im  8(/sc  —  Im) 

where 


(77) 


where  the  points  (V),  /,),  i = 1,2  are  obtained  from  the  I-V  char¬ 
acteristics,  considering  the  value  of  A /  so  that  the  parallel  drawn 
through  I 1  to  the  Ox  axis  intersects  the  characteristic  for  the 
highest  illumination  level  just  a  little  below  the  maximum  power 
point,  see  Fig.  4. 

Pysch  et  al.  proposed  a  modification  where  at  least  three  I-V 
characteristics  were  used  to  determine  the  series  resistance.  These 
I-V  characteristics  can  be  obtained  for  three  levels  of  illumination: 
one  sun,  slightly  above  and  slightly  below  one  sun  [40], 

Method  26.  Araujo  and  Sanchez  developed  the  area  method  to 
determine  the  series  resistance  of  a  solar  cell  [29].  In  this  method 
the  one  diode  model  was  used,  Eq.  (3).  Rsh  is  considered  infinite 
and  the  ideality  factor  of  the  diode  has  the  value  equal  to  one.  In 
this  case,  the  series  resistance  can  be  calculated  as  follows: 


Rs  =  2 


(Voc 

[isc 


A 

t 


(75) 


where  A  represents  the  area  between  the  I-V  characteristic  curve 
and  the  I  and  V  axes. 

Method  27.  Cape  et  al.  developed  the  flash  lamp  method  [41], 
This  method  allows  the  determination  of  the  series  resistance 
using  Eq.  (76).  The  one  diode  model  was  used  and  Rsh  is  considered 
infinite. 


Rs  =  Rl 


(76) 


D  _  Um/{Isc—Im)]+  In  [(/sc — 7m) / fsc] 

2Vm-V0C 


(78) 


Method  29.  The  method  developed  by  Cotfas  et  al.  allows  the 
determination  of  the  series  resistance  using  the  one  diode  model 
[44].  This  method  uses  two  I-V  characteristics  -  one  is  measured 
and  the  other  is  ideal,  Fig.  5.  For  the  ideal  case  the  I-V  character¬ 
istic  can  be  calculated  with  n=  1  or  with  the  real  values  for  n  and  I„ 
obtained  with  method  21.  The  series  resistance  can  be  calculated 
with  Eq.  (79). 

Rs=i0.=  Vl*°)-Vm  (79) 

lm  'm 


where  Videai  is  obtained  on  the  ideal  I-V  characteristic  for  Im. 

Method  30.  Agarwal  et  al.  developed  a  method  which  allows 
the  determination  of  the  series  resistance  using  the  one  diode 
model  and  the  variation  of  photogenerated  current  function  of  the 
levels  of  illumination  [45].  The  shunt  resistance  for  this  method  is 
considered  infinite.  The  short  circuit  current  can  be  approximated 
with  the  photogenerated  current  for  low  illumination  levels, 
smaller  than  750  W/m2.  Under  low  illumination  levels,  it  shows 
linear  dependency  vs.  the  illuminations  levels.  For  illumination 
levels  higher  than  750  W/m2  the  dependency  is  under  linear. 
The  difference  between  the  photogenerated  current,  by  projecting 
the  tangent  of  the  Isc  at  low  7r,  with  the  Isc  curve  obtained  at  normal 
/t,  becomes  significant.  The  series  resistance  can  be  found  using 
the  slope  of  the  semi-logarithmical  dependence  ln(Iph-Isc)  vs.  Isc 
which  is  linear  for  illumination  levels  above  750  W/m2: 


where  RL  represents  the  load  resistance  and  VL  denotes  the  voltage 
across  the  load  resistance.  This  method  uses  two  very  short  pulses 
with  duration  of  1  ms,  to  determine  Voc  and  VL.  The  open  circuit 
voltage  is  determined  by  using  the  first  light  pulse.  The  voltage 
across  the  load  resistance,  which  was  previously  measured  with 
precision,  is  determined  using  the  second  light  pulse. 

Method  28.  The  maximum  power  point  method  which  allows 
the  determination  of  the  series  resistance  was  developed  by 
Picciano  [42,43],  The  method  uses  the  one  diode  model  and  the 
I-V  characteristic  is  measured  for  a  single  illumination  level. 
The  shut  resistance  is  considered  infinite  and  the  ideality  factor 
of  the  diode  has  a  constant  value  over  the  entire  I-V  characteristic. 


Fig.  4.  The  method  with  two  I-V  characteristics  for  the  determination  of  solar  cell 
series  resistance. 


In  =R 

V  lo  J  nVT 


(80) 


where  n  and  /„  can  be  calculated  using  method  21  or  can  be 
approximated  with  values  from  ideal  cases. 

Method  31.  The  series  resistance  of  solar  cells  can  be  deter¬ 
mined  using  the  method  developed  by  Polman  et  al.  [46],  This 
method  uses  the  two  diode  model  for  the  solar  cell  description. 
The  shunt  resistance  and  the  ideality  factor  of  diode  are  taken  into 
consideration.  Using  Eq.  (81),  the  series  resistance  can  be  calcu¬ 
lated  as  the  average  between  the  values  obtained  when  n  takes  the 
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Fig.  5.  The  Cotfas  method. 
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values  1  and  2. 


where  (dI/dV)v  =  v  is  obtained  from  the  slope  of  the  linear  fating 
of  a  few  points  around  V=Voc. 

Method  32.  The  method  proposed  by  Cowley  and  Sze  to 
determine  the  series  resistance  consists  of  measuring  the  semi- 
logarithmical  I-V  characteristic  under  dark  conditions  and  using 
the  one  diode  model.  The  series  resistance  can  be  calculated  using 
the  gap  value  on  the  V  axis,  between  the  semi  logarithmical  I-V 
characteristic  and  the  diffusion  line  [39,47], 

Method  33.  The  method  allows  the  determination  of  the  series 
resistance  of  the  solar  cells  using  both  measurements  under  both 
dark  and  illumination  conditions.  Aberle  et  al.  [48]  proposed  the 
determination  of  the  series  resistance  using  Eq.  (82): 


Y  V _ V dark.m  ^  m 

Im  Im 


(82) 


where  Vuark,m  represents  the  voltage  measured  in  dark  conditions, 
corresponding  to  the  current  Im  measured  in  illumination  condi¬ 
tions.  Because  the  effect  of  dark  series  resistance  was  neglected, 
the  error  can  be  larger  than  5%.  The  series  resistance  is  calculated 
with  the  corrected  equation  [49]  in  order  to  minimize  the  errors: 

_  Vdark'in  —  Vm— Use— Im)Rs  dark  (83) 

Im 


where  RSidark  can  be  calculated  using  Eq.  (53). 

Method  34.  The  ideality  factor  of  the  diode  can  be  found  using 
the  analytical  method  developed  by  Bayhan  and  Kavasogl  [50],  In 
this  case  the  measurements  are  performed  under  dark  conditions; 
the  authors  used  the  one  diode  model  for  the  solar  cell  description 
and  the  Lambert  W-function.  The  value  of  n  is  calculated  for  two 
regions  of  I-V  characteristic  measured  under  dark  conditions.  The 
first  region  is  for  the  voltage  interval  (0.2,  0.6)  and  the  second 
region  is  for  the  voltage  interval  (0.6,  0.8). 


3.  Results  and  discussion 

This  article  outlines,  discusses  and  compares  a  number  of  34 
theoretical  and  experimental  approaches  in  order  to  determine 
the  important  dc  parameters  of  solar  cells,  like  Rs,  Rsh,  n,  Ia  etc.  The 
analysis  of  the  I-V  characteristic  is  an  important  tool  for  the 
determination  of  the  above  quantities  and  it  can  be  measured 
under  both  illumination  and  darkness  conditions.  Most  methods 
use  the  I-V  characteristic  obtained  under  illumination  conditions 
because  these  are  real  operating  conditions  for  a  solar  cell.  The 
methods  that  use  the  dark  measurements  are  few,  for  example 
methods  14, 15, 16,  32  and  34.  There  is,  also,  a  method  which  uses 
both  the  dark  conditions  and  the  illumination  conditions,  namely 
method  33  developed  by  Aberle  et  al. 

The  equation  which  mathematically  describes  the  behavior  of 
solar  cells  depends  on  the  model  chosen  for  the  solar  cell.  The 
most  widely  used  model  is  the  one  diode  model  because  it  is  the 
simplest  and  describes  sufficiently  well  the  characteristic  of  most 
solar  cells  behavior,  as  it  can  be  seen  for  all  methods  except  the 
methods  8, 11, 15, 16  and  31. 

The  first  sixteen  methods  presented  in  this  paper  allow  the 
determination  of  all  the  parameters:  Iph,  Ia,  Rs,  Rsh  and  n.  The  values 
of  he,  Voc,  Pm  and  implicitly  Vm  and  Im  can  be  determined  from  the 
experimental  data  [10].  Two  or  more  parameters  of  a  solar  cell  are 
determined  using  methods  17-24.  The  last  nine  methods  allow  for 
the  determination  of  one  parameter,  especially  the  series  resis¬ 
tance  which  is  the  mostly  tried  parameter  for  solar  cells. 


Most  of  the  methods  which  allow  the  determination  of  all  the 
parameters  use  the  least-squares  numerical  techniques.  Eqs.  (3)  and 
(4)  are  implicit  equations.  The  non-linear  fitting  procedure  becomes 
complicated.  Eq.  (3)  with  the  approximation  RSI<V  becomes  an 
explicit  equation  and  the  fitting  procedure  becomes  simpler.  The 
above  approximation  holds  for  a  solar  cell  with  small  area  which 
implies  that  the  current  is  small  or  for  low  level  illuminations.  Using 
the  Lambert  W-function,  Eq.  (3)  becomes  an  explicit  equation  I=I(V) 
[51]  and  the  fitting  procedure  simplifies  the  approach.  Moreover, 
using  the  assumption  given  by  Eq.  (15),  the  explicit  equation  I=I(V) 
has  only  three  fitting  parameters.  The  extracted  process  is  easier  using 
the  consecrated  software  such  as  Matlab,  Mathematica  and  others.  In 
Matlab  only  a  few  lines  of  code  are  necessary  [17]. 

The  accuracy  of  the  fitting  procedure  depends  on  the  following 
factors:  the  choice  of  the  fitting  algorithm,  especially  on  the  initial 
values  of  the  fitting  parameters  and  the  confidence  interval.  Devel¬ 
opment  of  the  generic  algorithm  and  the  particle  swarm  optimization 
may  improve  the  accuracy  of  the  solar  cell  parameter  extraction.  A 
high  advantage  of  the  method  14  in  comparison  with  the  fitting 
methods  is  that  even  without  very  good  choice  of  the  values  of  the 
initial  fitting  parameters,  the  final  results  are  very  good  and  the  errors 
are  very  small  in  comparison  to  other  experimental  errors  like  those 
of  methods  25,  26  etc.  In  method  10  the  initial  values  of  the  fitting 
parameters  are  firstly  calculated,  leading  to  improvements  of  the 
results  obtained  by  the  fitting  process. 

To  solve  the  statistical  instability  in  the  fitting  algorithm  and  /- 
V  characteristic  in  concentrator  solar  cells,  Araki  and  Yamaguchi 
proposed  a  new  mathematical  model  which  describes  the  solar 
cell  behavior: 


,  ,  ,  ,  ^  V+IRS 

- ^ - +e 


(84) 


where  e  is  an  error  term  distributed  through  a  Gaussian  pdf 
centered  at  zero  [52], 

The  area  method  and  the  generalized  area  method  introduce  a 
significant  uncertainty  brought  in  through  the  area  limited 
between  the  two  axes  and  the  I-V  characteristic.  Sometimes,  the 
result  for  Rs  comes  out  negative  [39], 

The  methods  which  use  few  approximations,  such  as  the 
methods  10,  13  and  14  have  to  be  applied  to  obtain  good  values 
for  the  solar  cell  parameters.  The  methods  in  which  the  ideality 
factor  of  diode  is  given  the  value  equal  with  1  and  the  shunt 
resistance  is  considered  infinite  must  be  used  carefully.  Method  29 
shows  how  the  value  of  the  series  resistance  changes  when  the 
value  of  n  takes  the  real  value  or  the  theoretical  value  (n  =  1). 

The  solar  cell  parameters  such  as:  Rs,  n,  Rsh  and  Ia  depend  on  the 
irradiance  levels  [17,40,44,53],  Therefore  the  methods  17,  21  and 
25  which  used  two  or  three  I-V  characteristics  obtained  for 
different  levels  of  illumination,  may  be  used  carefully.  The  two 
characteristics  method  25  can  be  improved  using  three  or  more 
I-V  characteristics  and  Rs  can  be  calculated  from  the  inverted  slope 
of  the  linear  dependency  of  the  (V,,  /,-),  i=  1,2,3,... 

The  series  resistance  obtained  in  dark  conditions  is  under¬ 
estimated.  The  series  resistance  is  higher  due  to  the  larger  lateral 
electron  flow  in  the  emitter  under  illuminated  conditions,  the 
normal  operating  mode  for  solar  cells. 


4.  Conclusions 

In  this  review  paper  34  methods  developed  over  the  past  35 
years  to  determine  the  main  electric  parameters  of  a  solar  cell 
were  critically  presented,  assessed  and  discussed.  The  main  para¬ 
meters  in  concern  are,  the  reverse  diode  saturation  current,  the 
ideality  factor  n  of  diode,  the  series  resistance,  Rs,  the  shunt 
resistance,  Rsh,  and  the  photocurrent,  lph.  The  methods  outlined 
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fall  into  three  main  groups.  The  first  group  includes  the  methods 
which  are  based  in  the  theoretical  analysis  of  the  l-V  curve  based 
on  techniques  which  range  from  simple  regression,  to  non-linear 
regression,  stochastic  process  and  fitting  based  on  genetic  func¬ 
tions.  This  family  of  approaches  may  provide  the  values  of  all  the 
above  basic  solar  cell  parameters  which  play  the  central  role  to  the 
solar  cell  performance  investigations.  The  second  group  includes 
methods  which  succeed  in  providing  the  values  of  most  of  the 
important  parameters,  as  the  ones  above,  by  dealing  with  graphi¬ 
cal  analysis  of  the  l-V  measured  in  various  illumination  conditions, 
either  at  transient  or  at  steady  state.  The  third  group  of  methods 
uses  both  sophisticated  theoretical  analysis  of  the  l-V  character¬ 
istic  based  on  the  one  or  two  diode  models  and  graphical  analysis 
of  the  l-V  curves. 

These  methods  are  deployed  in  such  a  way  that  they  may  be 
effectively  combined  in  the  most  efficient  way  as  to  provide  the 
values  of  the  above  five  electric  parameters  of  a  solar  cell.  As  result 
of  the  assessment  of  the  34  methods  outlined  in  this  paper,  the 
researcher  may  choose,  according  to  the  experimental  set  up  and 
the  tools  possessed,  some  of  the  inter-related  approaches  outlined 
above,  in  order  to  get  a  set  of  values  concerning  the  above  solar 
cell  parameters  useful  for  further  investigation. 
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